Phenomenology of Photoemission Lineshapes of High Tc Superconductors 
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We introduce a simple phenomenological form for the self-energy which allows us to extract important 
information from angle resolved photoemission data on the high Tc superconductor Bi2212. First, we find 
a rapid suppression of the single particle scattering rate below Tc for all doping levels. Second, we find 
that in the overdoped materials the gap A at all k-points on the Fermi surface has significant temperature 
dependence and vanishes near Tc- In contrast, in the underdoped samples such behavior is found only at 
k-points close to the diagonal. Near (7r,0), A is essentially T- independent in the underdoped samples. The 
fiUing-in of the pseudogap with increasing T is described by a broadening proportional to T — Tc, which is 
naturally explained by pairing correlations above Tc- 
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Angle-resolved photoemission spectroscopy (ARPES) 
has played a major role in developing our understanding 
of the high Tc superconductors. The momentum and fre- 
quency resolved information contained in the one-particle 
spectral function 0] probed by ARPES provides critical 
insights difficult to obtain from other techniques. How- 
ever, an important open problem in interpreting ARPES 
data is the absence of simple representations of the spec- 
tral lineshape, analogous to the Drude formula for optical 
studies. 

In this paper, we take a step in this direction by in- 
troducing a simple phenomenological form for the self- 
energy which captures much of the important low fre- 
quency information contained in the ARPES data for the 
normal, superconducting and pseudogap phases of high 
Tc superconductors. We first test this on ARPES data 
at \iF for overdoped Bi2212, from which we determine 
the T-dependence of the superconducting (SC) gap and 
the one-particle scattering rate, the latter found to de- 
crease rapidly in the SC state, leading to the appearence 
of sharp quasiparticles at low temperatures. 

We then turn to underdoped Bi2212, where a highly 
anisotropic pseudogap is known to persist well above Tc 
The magnitude of this gap and its smooth evo- 
lution through Tc establish a strong connection between 
the normal state pseudogap and the SC gap below Tc, 
and suggest that the pseudogap arises from pairing cor- 
relations in a state without long range phase coherence 
1^. Here, we use our simple self energy to gain insight 
into the ARPES lineshape in the pseudogap regime. We 
find that near (tt, 0) A is T-independent, and the pseu- 
dogap fills in due to a T-dependent broadening which is 
naturally explained by pairing correlations above Tc- 

Finally, we use our analysis to shed more light on the 
very recent surprising result that the pseudogap at differ- 
ent k-points turns off at different temperatures, leading 
to gapless arcs above Tc which expand in length until the 



entire Fermi surface is recovered at T* We show be- 
low that the T-dependence of the underdoped lineshape 
away from (tt, 0) is rather different, with the vanishing of 
the pseudogap controlled instead by a T-dependent A. 

The data analysed below have been reported earlier 
[|,|§. The ARPES intensity /(k,tj) is proportional 
to /(cj)yl(k, w), where / is the Fermi function and A 
the spectral function Recently, we have proposed 
a method |^ which allows us to eliminate the effect 
of / from ARPES data and focus directly on A. In 
brief, using the mild assumption of particle-hole sym- 
metry A{—e]i, ~uj) — A(ek,w) for small j^l and within 
the small k- window centered at k^?, one can show that 
the symmetrized intensity I{uj) + I{—uj) at k.p is sim- 
ply the spectral function (convolved with the resolution) 
. Results obtained from symmetrized data agree with 
those obtained from the leading edge of the raw data [|j . 
Here, we use the symmetrized data for two reasons: it is 
a useful visual aid, and (due to absence of Fermi cutoff) 
it allows more stringent comparisons to the fits. We have 
checked that the fits discussed below agree equally well 
with the raw data. 

We begin with the overdoped samples, where there are 
no strong pseudogap effects. The simplest self-energy 
which can describe the data at all T is 

E(k,a;) = -iTi + AV[(c^ + «0+) + e(k)]. (1) 

Here Fi is a single-particle scattering rate taken, for sim- 
plicity, to be an w-independent constant. It is effectively 
an average of the (actual w-dependent) E" over the fre- 
quency range of the fit. The second term is the BCS 
self-energy (corresponding to the diagonal term of the 
Nambu-Gorkov propagator) where A is the SC gap, and 
e(k) the dispersion (with e(kF) = 0). We emphasize that 
this is a minimal model; modeling the full w-dependence 
of the self-energy will involve more than Fi and A. 
In Fig. la, we show symmetrized data for an over- 
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doped Tc=82K sample at hp along (7r,0) — (tt, tt) 
together with the fits obtained as follows. Using 
Eq. 1 we calculate the spectral function 7rA(k, ut) — 
E"(k, to)/ [{u - ek - S'(k, + I]"(k, , convolve it 
with the experimental resolution, and fit to symmetrized 
data. The fit is restricted to a range of ±45 meV given 
the small gap in the overdoped case and the sharpness 
of the quasiparticle peaks below Tc- We find that Eq. 1 
describes the data quite well ||] . 

The T-variation of the fit parameters A and Fi are 
shown in Fig. 2a. A(T) decreases with T, and although 
small at Tc, it only vanishes above Tc, indicating the 
possibility of a weak pseudogap. This effect is sample 
dependent, in that several overdoped samples we have 
looked at, the gap vanishes closer to Tc- We caution that 
the error bars shown in Fig. 2a are based on the RMS er- 
ror of the fits, but do not take into account experimental 
errors in /i and k^?. 

ri(T) is found to be relatively T-independent in the 
normal state. In comparing this to optical conductiv- 
ity data, it must be kept in mind that Fi is an average 
of the true S" near (tt, 0) over the w-range of the fit, 
whereas 1 /rtr is a zone average weighted by velocity fac- 
tors which also differs from S" due to different thermal 
factors. Despite this, it is interesting to note that opti- 
cal conductivity work has also found a relative lack of T 
dependence to l/Ttr{uj) for uj > 2A Q. Below T^ we see 
that Fi decreases very rapidly, and can be perfectly fit to 
the form a + bT^ This rapid drop in linewidth lead- 
ing to sharp quasiparticle peaks at low T, which can be 
seen directly in the ARPES data, is consistent with SC 
state microwave and thermal conductivity pl] | measure- 
ments, and implies that electron-electron interactions are 
responsible for Fi [|^,^. Note the clear break in Fi at Tc, 
despite the fact A has not quite vanished. We have seen 
similar behavior to that described above for a variety of 
overdoped samples at several k points. 

We next turn to the more interesting underdoped case. 
We find that near (tt, 0) the self-energy (1) cannot give 
an adequate description of the data, in that it does not 
properly describe the pseudogap and its unusual "filling 
in" above T^. Theoretically, we cannot have a divergence 
in E(ki?,a; = 0) in a state without broken symmetry. A 
simple modification of the BCS self-energy rectifies both 
these problems: 



I](k, Lj) = -iVi + AV[cj + e(k) + iVo 



(2) 



The new term Fq (T) should be viewed as the inverse pair 
lifetime; below Tc, where the pairs have infinite lifetime, 
Fq = 0, and reduces to (^. The theoretical motiva- 
tion for (^) is given in the Appendix. We stress that this 
three parameter form is again a minimal representation 
of the pseudogap self-energy. Since it is not obviously a 
unique representation, it is very important to see what 
one learns from the fits. 



In Fig. lb, we show symmetrized data at the (tt, 0) — 
(tt, tt) Fermi crossing for a Tc=83K underdoped sample. 
Below Tc we see quasiparticle peaks. Above Tc these 
peaks disappear but there is still a large suppression of 
spectral weight around uj=0. As T is raised further, the 
pseudogap fills in (rather than closing) leading to a flat 
spectrum at a temperature of T* (200K). The self-energy 
gives a good fit to the data. These fits were done be- 
low Tc over a larger energy range (±75 meV) than in the 
overdoped case because of the larger SC gap. The range 
above Tc was increased to ±85 meV so as to properly 
describe the pseudogap depression. 

In Fig. 2b, we show the T-dependence of the fit pa- 
rameters. We find a number of surprises. First, A is 
independent of T within error bars. Similar behavior has 
been infered from specific heat and tunneling 
data. This T-independence is in total contrast to the be- 
havior of the overdoped 82K sample with almost identical 
Tc at the same k point. In addition, for the underdoped 
sample, the gap evolves smoothly through Tc- 

The single-particle scattering rate Fi(r) for the un- 
derdoped sample is found to be qualitatively similar 
to the overdoped case. It is consistent with being T- 
independent above Tc, but with a value over twice as 
large as the overdoped case (allowing Fi to vary above 
Tc does not improve the RMS error of the fits). Second, 
we see the same rapid decrease in Fi below Tc as in the 
overdoped case |lOj. Note again the clear break at Tc- 

The most interesting result is Fo(r). We find Fq = 
below Tc and proportional to T—Tc above. This behavior 
is robust, and is seen in all the fits that we have tried. 
Moreover, a non-zero Fq is needed above Tc to obtain a 
proper fit to the data (its effect cannot be reproduced 
by varying the other parameters). The fact that this T- 
dependence is exactly what one expects of an inverse pair 
lifetime (with a prefactor about twice the weak coupling 
BCS value; see Appendix), is a non-trivial check on the 
validity of the physics underlying Eq. 2. Further, we 
observe from Fig. 2 that T* corresponds to where A(T) ^ 
Fo(r). This condition can be understood from the small 
Lu expansion of Eq. 2. 

We note, in passing, that both ARPES [p^|jl5| and 
tunneling measurements indicate that the low tem- 
perature SC gap A(0) in Bi2212 increases as the doping 
is reduced. This increase closely tracks that of the T* 
values determined from ARPES |15|. This is further ev- 
idence linking T* to the onset of pairing. 

The next important question is whether the T- 
dependence near (tt, 0) described above exists at other 
k^r-points. To answer this, we have looked at T- 
dependent data for a number of underdoped samples at 
two different k vectors joj. All data at the (tt, 0) — (tt, tt) 
Fermi crossing give results similar to those for the 83K 
sample. However at the second k-point, about halfway 
between the first and the d-wave node along the (0,0) — 
(tt, tt) direction, we see quite different behavior. We 
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demonstrate this in Fig. 3a where symmetrized data for 
a 77K underdoped sample are shown. For the first k 
point, one clearly sees the gap fill in above Tc, with lit- 
tle evidence for any T-dependence of the position of the 
spectral feature defining the gap edge, just as for the 
83K sample. In contrast, at the second k point, the gap 
is clearly closing, indicating a strong T-dependence of 
A. Similar behavior is seen in other underdoped samples 
with Tc between 75 and 85K. 

In Fig. 3b, we show the T-dependence of A obtained 
from fits (over a range of ±66 meV) at the second k 
point for the 77K sample. A is found to be strongly T- 
dependent, being roughly constant below Tc, then drop- 
ping smoothly to zero above. The strong T-dependence 
of A makes it difficult to unambiguously determine Fq 
from the fits at this k-point. On theoretical grounds, 
we expect that, here too, there is a non-zero Fq, and 
the closing of the pseudogap is again determined by 
A(T) ~ Fo(T), however this condition is satified by the 
rapid drop in A(T), rather than the rise in Fo(T). For 
completeness, we also show A(T) for this sample at the 
(tt, 0) — (tt, tt) Fermi crossing, which has a similar behav- 
ior to that of the 83K sample. 

We see that these results give further evidence for the 
unusual k-dependences first noted in ref . . Strong pair- 
ing correlations are seen over a very wide T-range near 
(tt, 0), but these effects are less pronounced and persist 
over a smaller T-range as one moves closer to the diago- 
nal. This is clearly tied to the strong k-dependence of the 
effective interaction and the unusual (anomalously broad 
and non-dispersive) nature of electronic states near (tt, 0). 
Some of these features are captured in recent theoretical 
studies of the pseudogap p^^-p^ . 

To summarize, we have introduced a simple phe- 
nomenological self-energy expression which helps us to 
analyze ARPES data and gain insight into the T- 
dependences of important parameters like the gap and 
the one-particle scattering rate. Perhaps the most inter- 
esting new result concerns the modeling of the pseudogap 
data in the underdoped cuprates, where we found a new 
lifetime effect above Tc proportional to T — Tc- We ar- 
gue that this is naturally explained as an effect of pairing 
correlations above Tc on the one-particle spectral func- 
tion. In fact, this qualitative observation argues against 
non-pairing theories of the pseudogap, in that a term 
proportional to T — Tc does not naturally appear in such 
theories. A second important result concerns the differ- 
ences between the pseudogap behavior in different parts 
of k-space. Near (7r,0) we found a constant A(T) with 
"filling in" of spectral weight with increasing T. Away 
from this region, the pseudogap closes rather than fills in. 
These observations impose very important constraints for 
the microscopic theory of high Tc cuprates. 

Appendix: The simplest way to motivate the second 
term of Eq. 2 is as follows. The BCS E of Eq. 1 may be 
seen as arising from a bare Green's function Go dressed 



by a pair "susceptibility" of the form A^(5(q)(5(u;) rep- 
resenting static long range order. If the pairs have a 
finite inverse lifetime Fq, then 5{u>) is broadened into a 
Lorentzian, leading to Eq. 2. 

We next sketch a more formal derivation, valid in 
the regime of "small fluctuations" which gives 

further insight into the form of of Fo(T). The low- 
est order graph is Go dressed by the fluctuation pair 
propagator L: E(k, e„) = -T /(c^q)i(q, Wi.)Go(q - 
k,cj,, - e„). Here T-^(q,w^) = A^o(e + a\uj^\ + q^\l) 
with e ~ (T - Tc)/Tc, a = n/STc, the coherence 
length ^0 ~ vp/Tc, and iVo is the density of states. 
e„ and Wi, are Fermi and Bose frequencies respec- 
tively and J{dq) = J d^q/{2TT)^ . Evaluating us- 
ing standard contour integration, we obtain I](k, e„) = 
(T/A^o)/(dq)i^i(q)i^2(q). Here Ti(q) = l/[e + q'Q is 
sharply peaked in q with a scale ^^/e. Before ana- 
lytic continuation |ze„| > ttT, and thus the g- variation 
of T2(q) = l/[ien + z8Tc(e + q^^D/ir + e(q - k)] is on 
the much larger scale of ^(7^. Thus we make the approxi- 
mation S ~ F2(0)(T/iVo) /(i(q)Fi(q). Finally, doing the 
Bose sum in (|A(r,t)p) = T /(rfq)i(q, t^i/), we get 
the "fluctuating" gap (|A|2) = [T /N^) j{dq)Fi{q). Sub- 
stituting this in E, and setting ie„ — > w + iO^, we obtain 
the second term in (|) with Fo(T) = 8(T - Tc)/7r. 

In T) = dimensions, which may be relevant in the 
vicinity of the dispersionless (tt, 0)-point, we do not need 
to make the approximation of pulling ^2(0) out of the q- 
integration, which justifies the result even for arbitrarily 
small uj. In any case, we emphasize that above derivation 
is used here only to motivate the form of Eq. 2, which is 
then used for fits in a more general context. 

Acknovifledgements; We thank Yuri Vilk, Andrei 
Varlamov, Boldizar Janko, and Jim Allen for useful dis- 
cussions. This work was supported by the U. S. Dept. 
of Energy, Basic Energy Sciences, under contract W- 
31-109-ENG-38, the National Science Foundation DMR 
9624048, and DMR 91-20000 through the Science and 
Technology Center for Superconductivity. 



[1] 
[2] 

[3] 
[4] 
[5] 

[6] 
[7] 



M. Randeria et al, Phys. Rev. Lett. 74, 4951 (1995). 

D. S. Marshall et al, Phys. Rev. Lett. 76, 4841 (1996); 

A. G. Loeser et al, Science 273, 325 (1996). 

H. Ding et al, Nature 382, 51 (1996). 

H. Ding et al, Phys. Rev. Lett. 78, 2628 (1997). 

For a detailed review and references, see: M. Randeria, 



Varenna Lectures, 1997; cond-mat/9710223 . 

M. R. Norman et al., preprint (1997); |cond-mat/9710163 



Spectra analyzed below were normalized to the constant 
signal above the Fermi energy, which was then subtracted 
off. No other background subtractions were made since 
we are focussing on low energy data. 



3 



[8] 

[9] 
[10] 

[11] 



[12] 
[13] 
[14] 
[15] 

[16] 
[17] 

[18] 
[19] 

[20] 



The deviations beyond about 40 meV at low T are due 
to the strong oj-dependence of the SC state self-energy 
for uj > 2A; see M. R. Norman et al., Phys. Rev. Lett. 
79, 3506 (1997). 

A. V. Puchliov, D. N. Basov, and T. Timusk, J. Phys. 
Cond. Matter 8, 10049 (1996). 

The residual term a is due to the finite uj range of the fit. 

Reducing this range leads to a reduction in a. 

D. A. Bonn, P. Dosanjh, R. Liang, and W. N. Hardy, 

Phys. Rev. Lett. 68, 2390 (1992); K. Krishana, J. M. 

Harris, and N. P. Ong, Phys. Rev. Lett. 75, 3529 (1995). 

J. W. Loram et al, J. Supercond. 7, 243 (1994). 

Ch. Renner et al., submitted, Phys. Rev. Lett. 

J. M. Harris et al, Phys. Rev. B 54, R15665 (1996). 

H. Ding et al., SNS97 proceedings, to be published, J. 

Phys. Chem. Solids. 

N. Miyakawa et al., to be published, Phys. Rev. Lett. 
V. B. Geshkenbein, L. B. loffe, and A. I. Larkin, Phys. 
Rev. B 55, 3173 (1997). 

X. G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). 
J. R. Engelbrecht, A. Nazarenko, M. Ra nderia, E. 



Dagotto, preprint (1997); cond-mat/9705166. 
E. Abrahams, M. Redi and J. Woo, Phys. Rev. B 1, 208 
(1970); A. Schmid, Z. Physik 231, 324 (1970); B. R. 
Patton, PhD thesis, Cornell Univ. (1971). 



100 



80 



60 



OD82K 





. • (a) 


• 
• 


• • • • • 


• 








m m ^ 

CD 

A 




14 1 T I 




- + 1 M T 



200 



UD83K 




20 





50 100 150 200 50 100 150 200 
T (K) T(K) 

FIG. 2. A (open circles), Fi (solid circles), and Fo (solid 
squares) versus T at the (7r,0) — (■7r,7r) Fermi crossing for (a) 
a rc=82K overdoped sample and (b) a rc=83K underdoped 
sample. The dashed line marks Tc. The error bars for A are 
based on a 10% increase in the RMS error of the fits. 
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FIG. 1. Symmetrized data for (a) a Tc=82K over- 
doped sample and (b) a rc=83K underdoped sample at the 
(tt, 0) — (tt, tt) Fermi crossing at five temperatures, compared 
to the model fits described in the text. 



FIG. 3. (a) Symmetrized data for a rc=77K underdoped 
sample for three temperatures at (open circles) k_F point 1 in 
the zone inset, and at (open triangles) Icf point 2, compared 
to the model fits, (b) A(r) for these two k points (filled and 
open circles), with Tc marked by the dashed line. 
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